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ABSTRACT 

In order to study systematically the soft X-ray emission of Active Galactic 
Nuclei (AGNs) at medium to high redshifts, we have analyzed ROSATVSPC 
and HRI data of QSOs at 0.26 < z < 3.43 selected from the second deepest 
ROSAT PSPC survey carried out in 1991-1993 by McHardy et al. (1998). 
Our sample of 22 type 1 QSOs is nearly complete above a flux limit of 1.4 x 
10~^* erg cm^-^ s^^ in the 0.1 — 2 keV band. Of these, nine QSOs show long 
term (~ 2 yr) X-ray variability by a factor of 1.5 — 3.5. Significant excess 
absorption above the Galactic column is seen in three QSOs. The soft X-ray 
photon index of the QSOs ranges from 1.4 to 3.7. Three QSOs have steep 
soft X-ray spectra {Tx > 3.0), one of which is a narrow- line QSO - a high 
luminosity version of narrow-line Seyfert 1 galaxies. The average photon index 
(< Fx >) is2.40±0.09 (with a dispersion of 0.57) in the 0.1-2 keV band. The 
average QSO spectra in four redshift bins flatten from an average photon index 
of < Fx >^ 2.53 at 0.25 < z < 1 to < Fx 2 at 2 < z < 3.4. The flattening 
of the average photon index can be understood in terms of the redshift effect 
of the mean intrinsic QSO spectrum consisting of two components - a soft X- 
ray excess and a power-law component. We have also studied optical spectra 
of 12 of the 22 QSOs. 
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1 INTRODUCTION 

QSOs or quasars represent the high luminosity end of a class of objects known as active 
galactic nuclei (AGNs). They are the most luminous continuously emitting objects in the 
universe and emit over the entire range of electromagnetic waves. It is widely believed that 
high energy emission (X-rays and 7-rays) from an AGN originates in the innermost region 
of accretion disk around a super-massive black hole (SMBH) and much of the low energy 
emission e.g., infrared, optical, ultraviolet, is due to the reprocessing of high energy photons 
in a medium surrounding the accretion disk. While the atomic line emission in the optical, 
UV, and X-ray probe the circumnuclear medium surrounding the SMBH, the shape of the X- 
ray spectrum described by photon index (Fx) is a critical parameter to constrain competing 
models for the mechanisms of X-ray continuum emission. 

The increased sensitivity of the position sensitive proportional counter (PSPC) on-board 
ROSAT compared to earlier missions allowed a significant improvement in the study of 
the soft X-ray emission of AGNs. The X-ray spectral shape of quasars has been studied 
extensively using earlier missions such as HEAO-I, Einstein, EXOSAT, and Ginga (e.g., 
Mushotzky 1984; Wilkes & Elvis 1987; Canizares & White 1989; Comastri et al. 1992; 
Lawson et al. 1992; Williams et al. 1992). Some of these missions were not sensitive below 

2 keV. These earlier studies suggested that the X-ray emission of quasars is well described 
by a power- law with photon index Vx ~ 1.5 for radio- loud quasars and ~ 2.0 for radio 
quiet quasars. Large samples of AGNs have been studied by Walter & Fink (1993), Wang, 
Brinkmann & Bergeron (1996), and by Rush et al. (1996) using the ROSATVSPC However, 
most of the objects studied in the above papers are nearby and intrinsically bright AGNs. 
Also the AGN samples studied in the above papers are not complete and the derived results 
may be biased by selection effects. Laor et al. (1994, 1997) studied X-ray spectra of a 
complete sample of 23 quasars from the Bright Quasar Survey (BQS) with z < 0.4 and 
Nh < 1.9 X 10^° cm~^. They found an average photon index of 2.63±0.07 for their complete 
sample. Reeves & Turner (2000) studied a larger sample of 62 quasars with redshift in the 
range 0.06 - 4.3 and My < —23. They found an average photon index. Fx = 1-66 ± 0.04 
for 35 radio-loud quasars and Fx — 1-89 ± 0.05 for 27 radio-quiet quasars in their sample. 
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Table 1. Basic pcirameters of QSOs 



MJM^ No. ROSAT name HRI position^ Optical position^ Radio* 









a(2000) 


<5(2000) 


a(2000) 


5(2000) 






2 


RX 


J1334.7+3800 


13 34 41.87 


38 00 11.4 


13 


34 


41.82 


38 00 11.3 


0.26 


18.69 


3 


RX 


J1333.7+3803 


13 33 42.41 


38 03 35.3 


13 


33 


42.36 


38 03 36.3 


1.069 


18.60 


7 


RX 


J1334.3+3757 


13 34 17.52 


37 57 22.1 


13 


34 


17.52 


37 57 22.4 


1.14 


18.35 


10 


RX 


J1334.2+3759 


13 34 10.57 


37 59 56.4 


13 


34 


10.62 


37 59 56.3 


0.38 


19.55 


11 


RX 


J1333.5+3746 


13 33 32.32 


37 46 42.4 


13 


33 


32.01 


37 46 41.1 


0.826 


20.74 


13 


RX 


J1333.9+3759 


13 33 58.49 


37 59 39.3 


13 


33 


58.55 


37 59 38.2 


1.61 


21.10 


15 


RX 


J1334.7+3759 


13 34 42.81 


37 59 15.0 


13 


34 


42.77 


37 59 15.0 


1.14 


19.83 


17 


RX 


J1334. 1+3754 


13 34 00.96 


37 54 04.7 


13 


34 


01.03 


37 54 04.9 


1.64 


21.12 


18 


RX 


J1335.7+3751 


13 35 44.46 


37 51 43.3 


13 


35 


44.66 


37 51 40.8 


1.621 


20.25 


20 


RX 


J1335.5+3757 


13 35 30.26 


37 57 50.5 


13 


35 


30.30 


37 57 50.0 


1.39 


20.86 


21 


RX 


.11334.5+3748 


13 34 31.30 


37 48 31.2 


13 


34 


31.33 


37 48 31.4 


1.359 


22.27 


23 


RX 


J1333.7+3757 


13 33 44.34 


37 57 53.6 


13 


33 


44.27 


37 57 52.6 


0.97 


20.37 


24 


RX 


J1335.6+3757 


13 35:35.48 


37:57:46.3 


13 


35 


35.48 


37 57 46.2 


1.63 


20.21 


29 


RX 


J1335.7+3755 


13 35 42.38 


37 55 43.6 


13 


35 


42.51 


37 55 41.8 


1.90 


19.30 


30 


RX 


J1334.9+3757 


13 34 52.23 


37 57 44.7 


13 


34 


52.16 


37 57 44.8 


1.89 


20.47 


31 


RX 


J1333.9+3752 


13 33 55.94 


37 52 57.5 


13 


33 


55.81 


37 52 58.6 


2.14 


20.32 


37 


RX 


J1334.4+3746 


13 34 24.43 


37 46 15.9 


13 


34 


24.57 


37 46 15.2 


1.570 


20.11 


48 


RX 


J1335.5+3804 


13 35 30.00 


38 04 31.5 


13 


35 


29.69 


38 04 32.7 


0.692 


19.17 


56 


RX 


J1334.7+3757 


13 34 45.43 


37 57 21.5 


13 


34 


45.35 


37 57 22.8 


1.89 


20.36 


57 


RX 


,11333.6+3754 


13 33 35.69 


37 54 16.7 


13 


33 


35.62 


37 54 13.2 


1.525 


21.68 


61 


RX 


J1333.6+3749 






13 


33 


34.86 


+37 49 16.92 


3.43 


19.73 


63 


RX 


J1334.4+3806 


13:34:22.44 


38:06:21.5 


13 


34 


22.24 


38 06 20.1 


2.593 


21.40 



(r) 



(r) 



'^Identification number in the catalog of McHardy et al. (1998). 
^Derived from ROSAT YiKl observation of 1997. 
^McHardy et al. (1998) 

(r) indicates that the source is detected in the preliminary 20-cm radio map. 



Although the sample of Laor et al. (1997) is complete, it is confined to the nearby universe. 
The sample of Reeves & Turner (2000) includes some high redshift quasars but most of the 
quasars are nearby and the sample is not complete. 

We present the results of a detailed spectral and timing analysis of a nearly complete 
sample of 22 type 1 QSOs obtained from the ROSAT PSPC deep survey of McHardy et 
al. (1998). The basic parameters of the QSOs are hsted in Table 1. We will refer to all 
the objects as QSO, independent of their luminosity, but note that some have low X-ray 
luminosity (see Table 4). The aims are to extend the study of soft X-ray properties of QSOs 
at higher redshifts and to investigate whether QSOs at high redshifts have excess soft X-ray 
emission similar to that seen in the narrow-line Seyfert 1 galaxies. The outline of the paper 
is as follows. In §2 we describe the sample; in §3 we describe the X-ray observations and 
analysis. Section 4 deals with optical spectroscopy. In §5, we compare our results with other 
studies and discuss some of the implications. Finally, we conclude our study in §6. 

Throughout the paper, luminosities are calculated assuming isotropic emission, a Hubble 
constant of Hq = 75 km s~^ Mpc~^ and a deceleration parameter of go = unless otherwise 
specified. 
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Figure 1. MJM QSOs observed with the ROSAT HRI. The image was binned by 16 pixels (8 arcsec) and then smoothed by 
convolving with a Gaussian of u = 2 pixels. 



2 THE QSO SAMPLE 

The QSOs were selected from the second deepest ROSAT PSPC survey of McHardy et 
al. (1998) which covers a circular region of radius 15 arcmin of the sky towards the direction 
where the Galactic Nh is low (~ 7.9 x 10-*^^ cm"^). The ROSAT survey consists of two 
observations with PSPC and one observation with HRI (see below). There are 32 QSOs 
identified to a fiux limit of ~ 2 x 10^^^ erg cm^^ s^^ in the energy band of 0.5-2.0 keV (see 
Jones et al. 1997; McHardy et al. 1998). In order to be able to derive the shape of the soft X- 
ray spectrum of the individual QSOs from the single observation with the longest exposure 
and thus to avoid any effect due to the variation of soft X-ray emission over long time scales, 
we have selected 22 QSOs and excluded the faintest 10 QSOs which are not detected above 
5(7 level in any of the two PSPC observations. The soft X-ray selected sample of 22 QSOs 
(type 1) is nearly complete to a fiux limit of ~ 1.4 x 10^^"^ erg cm^^ s^^ in the energy band 
of 0.1—2.0 keV (corresponding to ~ 5 x 10^^^ erg cm^^ s^^ in the 0.5 — 2 kcV energy band). 
There are two X-ray sources, MJM 9 and MJM 54 brighter than the above flux limit but the 



Soft X-ray properties of a high redshift sample of QSOs observed with ROSAT 5 



— 1 1 r- 



— 1 1 r- 



10 - 



8 



m 
O 

CO 

cy 6 



° A 



_i I i_ 



_i I i_ 



Redshift (z) 



Figure 2. Redshift distribution of the 22 MJM QSOs 



optical nature of these objects is not known (McHardy et al. 1998) due to either lack of an 
optical counterpart or ambiguity in the identification. The ROSATYiBl image of the circular 
region of radius 15 arcmin is shown in Figure 1. The HRI image was binned by 16 pixels 
(8 arcsec) and then smoothed by convolving with a Gaussian of cr = 2 pixels. All the 22 
QSOs are marked in Fig. 1. The centres of the small circles are the optical positions of the 
QSOs as identified by McHardy et al. (1998). X-ray sources shown in Fig. 1, were classified 
as QSOs based on their optical spectra. All the 22 QSOs show broad (full width at half 
maximum (FWHM) > 1000 km s~^) permitted emission lines (Jones et al. 1997; McHardy 
et al. 1998). The hardness ratio (HR), defined as HR — where H and S are the count 
rates in the 0.5-2.0 keV and 0.1-0.5 keV bands, of the QSOs spans a range of ~ 0.14 — 1.78. 
In Figure 2, we show the distribution of the QSOs in the sample as a function of redshift 
and R mag. The number of QSOs increases with magnitude up to i? 20.5 mag and falls at 
fainter magnitudes (see Fig. 2). The QSO sample is optically complete at i? < 21 but it is 
likely that a few fainter QSOs have not been detected (McHardy et al. 1998). The redshift 
of the QSOs ranges from 0.26 to 3.43. 
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Table 2. Details of the ROSAT observations of the QSOs. 



Serial 


Sequence 


Instrument 


Start Time 


End Time 




Exposure 


No. 


No. 




Y, M, D, UT 


Y, M, D, UT 




Time (s) 


1. 


RH900717NOO 


HRI 


1997 06 04 16:12:58 


1997 07 13 22:26 


43 


201513 


2. 


RP900626N00 


PSPC 


1993 06 19 22:24:46 


1993 07 16 23:06 


28 


37658 


3. 


RP700283N00 


PSPC 


1991 06 23 20:59:40 


1991 06 26 20:44 


16 


71803 



3 X-RAY OBSERVATIONS AND ANALYSIS 

3.1 X-ray Observations 

The region of the sky containing the QSOs was observed twice with the ROSAT (Triimper 
1983) Position Sensitive Proportional Counter (PSPC) during 1991-1993, and once with 
the High Resolution Imager (HRI) (Pfeffermann et al. 1987) in 1997 June 19-July 16. The 
two PSPC observations carried out during 1991 June 23-26, and 1993 June-July together 
comprise the second deepest ROSAT PSPC survey. These observations were targeted at 
a(2000) = 13''34'"37.0", 5(2000) = +37°54'44" in the sky with an extremely low obscuration 
due to matter in our Galaxy (Nh = 7.9 x 10^^ cm~^). The PSPC deep survey data has 
already been reported by McHardy et al. (1998). The HRI observation was also a deep 
survey (exposure time = 201513 s) and was carried out in the same region of sky as the 
PSPC deep surveys. The details of the ROSAT observations are given in Table 2. All the 
QSOs lie within 15' from the field centre. 

3.2 Optical Identification 

The X-ray sources, presented here, were optically identified with QSOs using ROSAT PSPC 
and optical images, and optical spectra by McHardy et al. (1998). Here we have confirmed 
their identification using high resolution HRI X-ray images (central full width half maximum 
of ~ 4") . The analysis of HRI X-ray images has been carried out using the PROS^ software 
package. X-ray images were extracted from the observations listed in Table 2 and were 
smoothed by convolving with a Gaussian of a— 2". We have overlaid the contours of the 
smoothed HRI image of each QSO onto the corresponding optical image. The optical R band 
images the QSOs were obtained from the 2.5-m Isaac Newton Telescope (INT) using the 
wide-field camera on the night of 1999 April 19. The pixel scale is 0.333 arcscc per pixel and 
the seeing (FWHM) was ~ 1.3". The overlays of X-ray contours onto the R-band images are 

^ The PROS software package provided by the ROSAT Science Data Center at Smithsonian Astrophysical Observatory. 
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shown in Figure 3.2. The QSO MJM 61 was not detected in the HRI observation, therefore 
it is not shown in Fig. 3.2. The overlays reveal that X-ray emission from some of the QSOs 
is contaminated with that due to neighbouring sources. In particular, the X-ray emission 
from the QSOs MJM 18, 20, 23, and 24 may be contaminated significantly, however, the 
dominant X-ray emission is that due to QSOs. In the following, we have attributed the X-ray 
emission seen with the PSPC to the individual QSOs due to the difficulty in separating out 
the small contribution of neighbouring sources in the lower resolution (FWHM ~ 25 arcsec) 
PSPC images. 

The HEASOFT vS.O.l package was used to analyze the PSPC X-ray data of all QSOs. 
We have estimated the X-ray count rates for all the QSOs in our sample from the PSPC 
observations of 1991 and 1993. The total source counts for each QSO were obtained in 
the energy band of 0.1-2.0 keV from the unsmoothed PSPC images using a circular region 
centered on the peak positions, and after subtracting the background estimated from 3-4 
nearby circular regions. Note that for a few objects e.g., MJM 2, MJM 15 the centres of 
the circular region were shifted by a few arcsec in order to avoid the contamination from 
the neighboring objects. The radii of circles, used to extract the source counts, vary from 
45 — 90 arcsec. This was necessitated due to the different off-axis position of individual 
QSOs and to avoid contamination with the nearby X-ray sources. The background counts 
were estimated separately for individual objects using 3-4 nearby circular regions of radii 
1.5 — 2.5 arcmin in order to avoid any spatial variation of the background intensity. Net 
source counts were derived by subtracting appropriately scaled background counts. The 
count rates of all the QSOs, derived from the observations of 1991 and 1993, are listed in 
Table 3. Also hsted in Table 3 are the hardness ratio (HR) values where HR is defined as 
the ratio of the source counts in the energy band of 0.5-2.0 keV to that in the energy band 
0.1-0.5 keV. The QSOs MJM 31, 56, and 63, are not detected during the observation of 
1993. 

3.3 X-ray Light Curves 

In order to investigate the time variability of soft X-ray emission from the QSOs, we have 
extracted light curves from the ROSAT PSPC observations. The light curves for the source 
and the background were extracted using the 'xselect' package with time bins of 2000 s and 
in the PSPC energy band of 0.1-2.0 keV containing all the X-ray photons falling within 
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Table 3. Observed count rates and hardness ratios of MJM QSOs 



MJM No.i 




1991 




1993 


Variability 




count rate 


HR 


count rate 


HR 


factor 




10-2 cnt 




10-3 cnt s-i 






2 




— 0.213± 


1 S 4 -1- 
-LO . 1 HI u. 00 


_n 1 07 -L 04R 

KJ.i-^t U.U1U 


^ 2 




_LU . ^ _L- U . / U 


—0 ^fizi -1- n Ci'^fi 

\J.O\J^ zc u.uoo 


Q no -L 1 nc: 


_n c:qi -L n nocr 
U.OO-L m u.uoo 








—0 ^9'=1 -1- n 0^9 


1 4 zlQ -1- S4 


— fi/17 -1- O/IQ 
u.u'ii m u.u'iy 


^ 3 


10 


-Ly.-Lu u.uu 


—0 7fiS -4- n 09^ 
u. ( uo zn u.u^ij 


1 Ql -L 


— 7Q4 -1- 040 
u. / yH: ZJZ U.U1U 








_n 91 fi -1- n 077 

U.^-LU U.U( / 


c; Qi -1- R4 


— ^^^4 -t- 1 1 4 








— u.^Az in u.uoo 


7 70 -t- Q'^ 
/ . ( u in u.yo 


9RQ -t- 1 Cl^ 




1 ^ 




fi97 -t- 0*^0 

— Kj.v^i in u.uoy 


K OR -t- fiR 
o.uo in u.oo 


fi91 -t- 1 97 


2 




S 27 + 41 


—0 244 ±0114 


fi OS ± n 73 

U.UO _1_ L». 1 tJ 


—0 247 ±0111 


~ 2 


18 


Q 1 s -U n Afi 

O.J-O ^ u.tu 




0.00 ^ U. ( 


—0 OR'i -1- n 1 S7 




20(*) 


3.95 ±0.41 


-0.485 ±0.087 


2.78 ±0.56 


-0.221 ±0.195 


- 


21 


1.66 ±0.37 


0.109 ±0.216 


3.16 ±0.61 


0.144 ±0.204 


~ 1.5 


23(*) 


9.25 ±0.63 


-0.675 ±0.054 


6.50 ±0.77 


-0.600 ±0.103 




24(*) 


3.56 ±0.41 


-0.331 ±0.127 


4.35 ±0.59 


-0.264 ±0.132 




29(*) 


6.97 ±0.63 


-0.632 ±0.085 


6.94 ±0.74 


-0.602 ± 0.095 




30(*) 


1.80 ± 0.36 


0.047 ± 0.156 


3.67 ± 0.58 


-0.269 ±0.154 


~ 1.8 


31 


1.58 ± 0.34 


0.251 ± 0.221 


< 1.7 






37 


3.01 ±0.45 


-0.687 ±0.125 


2.83 ± 0.76 


0.638 ± 0.244 




48 


4.95 ± 0.51 


-0.532 ± 0.095 


2.54 ± 0.83 


-0.021 ± 0.325 


~ 2 


56 


2.83 ± 0.41 


-0.600 ±0.120 


< 0.65 






57(*) 


4.03 ±0.47 


-0.617 ±0.108 


3.49 ± 0.77 


-0.305 ± 0.203 




61 


4.23 ±0.58 


-0.372 ±0.135 


3.16 ±0.75 


-0.381 ± 0.213 




63(*) 


2.86 ± 0.40 


-0.436 ±0.146 


< 1.56 







^ (*) indicates that the X-ray emission form the source is contaminated by that due to neighbouring sources. The 
contamination could be as high as ~ 50% in case of MJM 20, 30, 55, 63, and 75. In other sources, the contamination is ~ 10% 
or less. 

"good time intervals". The source regions and the background regions were the same as 
described above. It was found that the background was variable by a factor of ~ 2 — 3 
during the observation of 1993 June 19. Therefore, the light curves obtained from the 1993 
observation are not suitable for variability studies. During the observation of 1991 June 
23, the background was reasonably constant. The background subtractions were carried out 
after appropriately scaling the background light curve to have the same area as the source 
extraction area. In order to test the constancy of the background light curves and soft X-ray 
variability of the QSOs, we have fitted constant count rates to the background light curves 
and source light curves. From the best-fit minimum values, it was found that background 
count rates obtained separately for each QSO are reasonably constant. It was also found 
that all QSOs, except MJM 10, do not show significant short-term variability during the 
observation of 1993 June-July. The soft X-ray variability characteristics of the QSO MJM 
10 has been discussed by Dewangan et al. (2001). 

A comparison of count rates during the observations 1991 and 1993 (see Table 3) reveals 
that 9 out of 22 QSOs show long term (time scale ~ 2 yr) variability. The QSOs MJM 2, 
15, 17, and 48 varied by a factor of 2 within 2 yr while the QSOs MJM 7, 21, and, 30 varied 
by factor of 3, - 1.5, and - 1.8, respectively The two QSOs MJM 56, and MJM 63, 
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Figure 3. Optical R band images of the QSOs overlaid by ROSAT HRI contours. In each optical image the central ob- 
ject is the QSO. The X-ray contours are drawn at intensity levels 0.119, 0.137, 0.155, 0.191, 0.245, 0.335, 0.515, 0.695, and 
0.965 count pixel"^ 

which are not detected above 3(7 level during the observation of 1993, varied by factor of 
more than 3.6 and 1.3, respectively. 

3.4 X-ray Spectral Analysis 

Photon energy spectra of all the QSOs in our sample were accumulated from their PSPC 

observation of 1991. The same regions for the source and the background, as described above 
(§4.1), were used. To improve the statistics, the ROSAT V'^VC pulse height data obtained 
in 256 pulse height channels for each QSO were appropriately re-grouped to have at least 
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20 counts in each channel. As examples, X-ray spectra of two QSOs thus obtained from the 
observation of 1991 are shown in Figure 4^. 

We used the XSPEC (Version 11.0) spectral analysis package to fit the data with spectral 
models. An appropriate response matrix and an auxiliary response file of the effective area 
of the telescope were used to define the energy response of the PSPC. The ROSAT PSPC 
spectra of QSOs shown in Fig. 4, were used for fitting spectral models. The X-ray spectrum 
of each QSO in the sample was first fitted with a redshifted power-law model with photon 
index, Fx, and absorption due to an intervening medium (Model A) with the absorption 
cross-sections as given by Balucinska- Church & McCammon (1992) and using the method 
of x^-minimization. The results of this fitting and the best-fit spectral model parameters are 
shown in Table 4. The errors quoted were calculated at the 90% confidence level based on 
X^njn-|-2.71. This simple model is a good fit to all the individual QSO spectra as evidenced 
by the minimum reduced values. The absorbing column densities derived by fitting the 
power-law model to the spectra of MJM 2, MJM 3, MJM 11, MJM 13, MJM 17, MJM 18, 
MJM 20, MJM 21, MJM 23, MJM 24, MJM 31, MJM 37, MJM 56, MJM 57, and MJM 61 are 
similar within errors to that due to our own Galaxy measured from 21-cm radio observations 
along the fine of sight to the source (Dickey & Lockman 1990), however, the error bars in 
the derived Nh are quite large for some objects. To better constrain the photon indices for 
these QSOs, we next fixed the absorbing column density to the Galactic value and fitted 
the power-law model (Model B) to each QSO. The results of this fitting are also given in 
Table 4. 

The absorbing column densities, 'iQ.lt^if x 10^*^ cm'^, 18.11^;? x 10^^ cm'^, 18.7l?:l x 
10^9 cm-^ derived for the QSOs MJM 7, MJM 10, and MJM 15 by fitting absorbed power- 
law model are in excess of the Galactic value (7.9 x lO^^cm"^) measured along their line of 
sight. This indicates that all the X-ray absorption may not be only due to matter in our 
own Galaxy but also due to matter local to the source. In order to verify the presence of 
excess absorption, we fixed the neutral hydrogen column density to the Galactic value and 
fitted the absorbed power- law model (Model B). In each case, although the fit is acceptable, 
there is an increase in the value of (A^^ = 5.33 for 23 dof in case of MJM 7; Ax^ = 
16.46 for 25 dof in case of MJM 10; Ax^ = 6.648 for 15 dof in case of MJM 15). To 
estimate the excess absorption local to each of the three QSOs, we introduced an additional 



^ Rest of the X-ray spectra can be made available on request to G. C. Dewangan (email: gulabQtifr.res.in) 



Soft X-ray properties of a high redshift sample of QSOs observed with ROSAT 11 



E 

o ^ 



Xo 

5 



MJM 2 (RX Jl 334.7+3800) 
• 1991 




MJM 3 (RX J 1333.7+3803) 
• 1991 




Energy (keV) 



0.5 

Energy (keV) 



Figure 4. Examples of the ROSAT PSPC spectra of MJM QSOs observed on 1991 and the best-fit models. The fitted model 
to each QSO spectrum is a redshifted power-law absorbed by our own Galaxy. 

component for absorption local to the source in our model B (Model C - redshifted power- 
law modified with Galactic as well as intrinsic absorption local to the source) and fitted to 
PSPC spectrum of the each of the three QSOs. The derived excess absorption, ANh, are 
9.71^^1° X lO^o cm-2 for MJM 7, 2.8tli x 10^° cm"^ for MJM 10, and 9.5^3.7 x lO^o cm'^ 
for MJM 15. The best-fit photon indices (Tx) are S.ltH, 3.7l|]i and 3.3l|];^ for MJM 7, 
MJM 10, and MJM 15, respectively. We measure the statistical significance of the reduction 
in the best-fit with addition of an intrinsic absorption component, described by Nh as a 
free parameter, using the F-test (Bevington 1969). The calculated values of F-statistic values 
and the corresponding probabilities are given in Table 5 for the QSOs MJM 7, 10, and 15. 
The F-statistics probability is the probability that the reduction in is not statistically 
significant. As can be seen the Table 5, at 95% level, the presence of additional absorption 
component local to the source is significant for the QSOs MJM 7, MJM 10 and MJM 15. 

For the QSOs MJM 29, MJM 30, MJM 48, MJM 61, and MJM 63, the amount of 
the absorbing column derived from the power-law fitting is unphysically small compared to 
the Galactic values of absorbing columns along their respective hues of sight. In order to 
constrain the photon indices, we fixed the absorbing column density to the Galactic value 
along the respective sight lines of individual QSOs and carried out the power-law model 
fitting to the PSPC spectra of the above QSOs. The PSPC spectra of the QSOs are well 
described by a power-law modified with absorbing column in our own Galaxy, as evidenced 
by the reduced xLin values (see Table 4) . 

In order to investigate any change in the soft X-ray spectral shape and flux, we have 
also analyzed the PSPC spectra of the QSOs observed in 1993. The PSPC spectra of the 
QSOs and the corresponding background spectra were extracted from the observation of 
1993 using the similar extraction regions as for the observation of 1991. The same spectral 
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Table 4. Spectral model parameters of QSOs derived from the observation 1991 
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"Model A(B) a Is redshlfted simple power-law model modified by an Intervening medium (Galaxy) described by Nh (Galactic 
Nj] ) at z=0.0. Model C is same as model B with an addition absorbing component at the source redshlft. 

''Observed flux in units of 10~i'l erg cm^2 §—1 in the energy band 0.1 — 2.0 keV. 

■^Unabsorbed soft X-ray luminosity in units of 10** erg s~i in the energy band 0.1 — 2.0 keV. 
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Table 5. F-test for the presence of excess absorption at z = zqso over the Galactic absorption of N/f = 7.9a;10^^cm ^ at 
z = derived for the 1991 observation 
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Table 6. Spectral model parameters of QSOs derived from the observation of 1993 
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"Models B and C are same as that defined in Table 4. 

''Observed flux in units of 10"^* erg cm-2 g-i in the energy band 0.1 — 2.0 keV. 
•^Unabsorbed soft X-ray luminosity in units of 10** erg s"^ in the energy band 0.1 — 2.0 keV. 



models, namely redshifted and absorbed power-law model with Nh as a free parameter and 
Nh fixed to the Galactic, as described above were fitted to the PSPC spectra obtained in 
1993. The power-law model with Nh as a free parameter is a good description to the observed 
PSPC spectra of individual QSOs. Except for the QSOs MJM 7 and MJM 10, the Galactic 
Nh along the line of sight of each QSO was found to be well within the range of derived 
Nh . Therefore, the power-law model with Nh fixed to the Galactic value was fitted to the 
spectrum of each QSO. The best-fit parameters derived are given in Table 6. 

The PSPC spectrum of the QSO MJM 7 is well described by a power-law of P = SMtHl, 
and Nh = 3.24;|:{^ x 10^° cm'^ (reduced xLn = 1-262 for 17 dof). The best-fit value 
of Nh is much higher than the Galactic value [N^""^ = 7.9 x 10^^ cm^^) suggesting excess 
absorption. The power-law fit with Nh fixed to the Galactic value is not acceptable (reduced 
Xmin — 2.123 for 18 dof). Addition of an absorbing component, apart from the Galactic 
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absorption component, local to the QSO gives an acceptable fit (reduced Xmin — 1-344 for 
17 dof) The best-fit parameters are Fx = 4.74ti;^?, and ANr = 2.861?;^? x lO^i cm'^. The 
F-statistic value for the addition of the intrinsic absorption component to the power-law 
with Nh fixed at the Galactic value is 11.4 and the F-statistic probability is 3.55 x 10~^. 
Thus the addition of an excess absorption at the QSO redshift to the power-law model with 
fixed Nh is a significant improvement at 99.4% level. 

The power-law model with Nh as a free parameter is an acceptable fit (reduced Xmin ~ 
0.83 for 11 dof) to the PSPC spectrum of MJM 10. The best-fit parameters are Nh = 
le.Olgi X 10^^ cm"2 and Fx = 3.7to:4- The best-fit value of Nr is significantly higher than 
the Galactic value (A^^"' = 7.9 x 10^^ cm~^) indicating excess absorption. The best-fit power- 
law model with Nh fixed to the Galactic value is also an acceptable fit (reduced xLin = 1-11 
for 12 dof) although the fit is becomes poor again indicating excess absorption. Next, we 
introduced an additional absorbing component local to the source and carried out the model 
fitting. The fit is an acceptable fit (reduced Xmin — 0-82 for 11 dof). The F-statistic value 
for the addition of new component to the power-law model with fixed Nh is 5.24 and the 
corresponding F-statistic probabihty is 4.28 x 10~^. Hence the power-law model modified 
by Galactic absorption and intrinsic absorption is a significant improvement at 90% level 
but not a significant improvement at 95% level over the power-law model with Galactic 
absorption. In the case of MJM 15, an excess absorption is not inferred from the 1993 
observations. This discrepancy could be due to the low S/N and poor spectral resolution of 
the PSPC spectra of MJM 15. 

3.5 Optical Spectroscopy 

Optical spectra of 12 of the 22 QSOs were obtained on the night of 1994 April 7 with the 
Multi Object Spectrograph (MOS) on the 3.6-m Canada France Hawaii Telescope (CFHT) 
as a part of a program to identify X-ray sources (McHardy et al. 1998). A 300 1 mm~^ 
grism in the first order with a Lorel3 CCD detector was used to cover a wavelength range 
of 4000 A-9000 A with ~ 15 A resolution. 

Optical spectra of the 12 QSOs were reduced using IRAF (v2.11.3). For details of the 
reduction see McHardy et al (1998). The fiux calibration is somewhat uncertain because the 
slitlets were not aligned at the parallactic angle. The fiux calibrated spectra were corrected 
for Galactic extinction by adopting the extinction law in Cardelli et al. (1989). Color excess 



Soft X-ray properties of a high redshift sample of QSOs observed with ROSAT 15 

Eb-v due to Galactic reddening along the line of sight of each QSO was calculated from 
the neutral hydrogen column density Nh using the relation 

Nh = 4.93 X 10^^ X Eb-v cm"^ (1) 

(Diplas & Savage 1994). The dereddened spectra were then deredshifted by dividing the 
starting wavelength and the pixel width in Angstrom by (1+z). The final spectra of 12 QSOs 
are shown in Figure 5. For each of the 12 objects, relative flux against the rest wavelength 
in A has been plotted. Strong emission lines have been marked in the figure. 

The analysis of the emission lines of MJM 10 has been discussed by Dewangan et 
al. (2001). Emission lines of other QSOs were modeled by Gaussian profiles using the IRAF 
task ngaussfit. As a first step, we fitted a single Gaussian profile to the strong emission lines. 
Most of the emission lines are described by a single Gaussian profiles, however, the broad 
wings of the Mg II hues in the spectra of the QSOs MJM 3, 7, and 20 are not well fitted by 
a single Gaussian profile. Therefore, the Mg II line in the spectra of the above QSOs were 
modeled by two Gaussian which gave satisfactory fit to the line core as well as to the wings. 
The two component fitting to the Mg II line revealed that the broad component is blue 
shifted with respect to the narrow component (see Table 7). The broad components of the 
Mg II line is blue shifted by ~ 1180, ~ 790, and ~ 1080 km s"^ in the spectra of MJM 3, 7, 
and 20, respectively. In the rest of the QSO spectra only one component is detected. This 
could be due to the low S/N ratios of the optical spectra. 

4 DISCUSSION 

4.1 The soft X-ray spectral shape 

The deep ROSAT PSPC observations together with improved sensitivity and energy resolu- 
tion, compared with earlier instruments, allows us to determine the shape of the soft X-ray 
continua of QSOs of redshifts up to 3.43. There are three QSO namely MJM 7, MJM 10, and 
MJM 15 which show steeper soft X-ray continua (F^ > 3.0). The QSO MJM 10 is known 
to be a narrow-line QSO (NLQSO) - high luminosity version of narrow-hne Seyfert 1 galax- 
ies (NLSls) (Dewangan et al. 2001). The average photon index of the QSOs is 2.40 ± 0.09 
with a dispersion of 0.57. 10 QSOs with z > 1.6 have mean photon index of 2.12 ± 0.10 
(with dispersion about mean of 0.39) and 12 QSO with z < 1.6 have mean photon index of 
2.62 ± 0.15 with a dispersion of 0.61. The average photon index of the QSOs with z > 1.6 
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Figure 5. Optical spectra of 12 of the 22 MJM QSOs. Vertical scale is the relative flux and the horizontal axis is the rest 
wavelength. The absorption feature seen in all the spectra is due to atmospheric absorption. 



in our sample appears to be flatter than the average photon index of QSOs at low redshifts. 
Walter & Fink (1993) analyzed the PSPC spectra of 58 AGNs observed in the ROSAT 
All Sky Survey (RASS) which contains a subsample of 24 QSOs with mean photon index, 
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Table 7. Emission line parameters of the QSOs. The quantities in brackets represent la errors. 
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Corrected for the instrumental broadening (FWHM~ 15 A). 



Fx — 2.57 ±0.06 (Laor et al. 1997). A similar average photon index of 2.65 ±0.07 was found 

by Schartel et al. (1996) for 72 QSOs from the LBQS sample detected in the RASS. Laor et 
al. (1997) studied the ROSAT PSPC spectra of a complete sample of 23 optically selected 
QSOs with z < OA and found average photon index of 2.62 ± 0.09 for their sample. 
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Figure 6. Best-fitting power-law plioton indices obtained by fitting simple power-law models to the average spectra of MJM 
QSOs at four redshift bins. Also plotted is the average photon index of 23 QSOs from the BQS sample of Laor et al. (1997). 
The average photon index of QSOs flattens at higher redshift. 



In order to further explore the redshift dependence of soft X-ray spectral slope of QSOs, 
we formed average spectra of QSOs in four redshift bins - z — 0.25 — 1.0 (5 QSOs), z — 
1.0 - 1.6 (7 QSOs), z = 1.6 - 2.0 (7 QSOs), and ^ = 2.0 - 3.43 (3 QSOs). The average 
spectra were then fitted by simple power-law modified by the Galactic absorption only. 
In Figure 6, we have plotted the average photon index as a function of redshift for the 
QSOs in our sample. Also plotted in Fig. 6 is the average photon index of 23 Quasars from 
the Bright Quasar Survey (BQS) sample of Laor et al. (1997). The average photon index 
flattens at higher redshifts e.g., < F^ >= 2.62 ± 0.09 at z < 0.4 and < F^ >= 2.17lo:o8 at 
1.6 < z < 2.0. 

Blair et al. (2000) studied soft X-ray spectral evolution with redshift using a large sample 
of 165 QSOs observed with RO SAT and found the 0.1 — 2 keV average spectra in 5 redshift 
bins harden from Fx ~ 2.6 at z — 0.4 to Fx ~ 2.1 at z — 2.4. This result is similar to 
that found here. Blair et al. (2000) also noted that the spectra in 0.5 — 2 keV band show 
no significant variation in spectral index with redshift suggesting the presence of a spectral 
upturn below 0.5 keV. While Blair et al. (2000) found that at lower redshifts spectra need a 
soft excess component in addition to a power-law and inclusion of a blackbody component 
{kT ~ 100 eV) for the soft excess results in a significant improvement over a single power- 
law, Laor et al. (1997) did not require such a soft excess component. The discrepancy could 
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be due to the fact that Laor et al. (1997) added 1% systematic errors in quadrature to the 
statistical errors. The difference between the average spectral slopes of QSOs at low and 
high redshifts can be understood in terms of mean X-ray spectrum consisting of a power- 
law of Vx ~ 1-9 and a soft excess below ~ 0.5 keV in the rest frame. At higher redshifts 
{z < 2) the rest frame soft excess component below ~ 0.5 keV is not covered in the observed 
0.1 — 2 keV band, the observed spectrum is a power-law with Vx ~ 1-9. While at lower 
redshifts, the observed spectrum, consisting of both soft excess and power-law component, 
is steeper when fitted with a single power-law and similar to that of z < 2 QSOs when the 
soft excess component below ~ 0.5 keV is excluded from the fit. This could also be the 
reason why Reeves & Turner (2000) found mean photon index of ~ 1.9 in the 0.6 — 10 keV 
ASCA band for their sample of 27 radio quiet quasars. 

4.2 Intrinsic Absorption 

Due to the low signal-to-noise of the PSPC spectra of the QSOs in our sample, we are not 
able to determine the consistency of HI columns derived from the 21 cm radio observations 
(Dickey & Lockman 1990) and HI columns derived from the absorbed power-law model fits 
to the observed PSPC spectra. Laor et al. (1997) have shown that the two columns agree 
to a level of about 5% — 9% at high Galactic latitudes. There are indications that some 
broad-line AGNs show intrinsic absorption (e.g.. Page et al. 2001; Akiyama et al. 2000). The 
power-law models modified by the Galactic absorption provide acceptable fits to the PSPC 
spectra of all QSOs except MJM 7, MJM 10, and MJM 15 (see Tables 4,& 5). Thus most of 
the QSOs in our sample lack intrinsic absorption. Only three QSOs, MJM 7, MJM 10, and 
MJM 15, out of 22 show indications for the presence of intrinsic absorption. The amount of 
intrinsic absorption inferred from the best-fit absorbed power-law models, ~ 10^^ cm~^ for 
MJM 7, is quite high and may cause the QSO to be redder than others in the sample. The 
ratio of soft X-ray to optical R band flux can be written as 

lo9{^) = log{fx) + ^ + 3.5 (2) 
fR 2.5 

where fx and fn are X-ray and optical R band fluxes, and tur is the R band magnitude. For 
MJM 7, the ratio of observed X-ray and optical flux (^) is calculated to be 2.95 x 10~^ while 
the ratio is 0.0160 for MJM 15 which is located at the same redshift (z=1.14) as MJM 7. In 
order to investigate whether the lack of soft X-ray emission of MJM 7 is only due to intrinsic 
absorption or the QSO is intrinsically X-ray weak, we calculated the ratio of unabsorbed 
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soft X-ray flux and observed R band flux for MJM 7 and MJM 15. Tlie ratio (^) is 0.017 
for MJM 7 and 0.072 for MJM 15. Tlius, tlie low observed X-ray flux of MJM 7 may not be 
entirely due to absorption and MJM 7 could be intrinsically weaker. The amount of intrinsic 
HI column in MJM 10 is small (2.8lo!8 ^ 10^° cm~^) and is similar to that observed in some 
NLSl galaxies (Grupe et al. 1998). 

4.3 Optical Spectra 

Optical spectra of the 12 QSOs, except MJM 10, are typical of QSOs. The optical spectra 
of MJM 10 is typical of NLSl galaxies or NLQSOs and has been discussed in detail by 
Dewangan et al. (2001). The other 11 QSOs show broad permitted lines in the ultra-violet 
region. The Mg II lines in the spectra of the QSOs MJM 3, 7, and 20 have been decomposed 
into a narrow and a broad component. The broad component is found to be blueshifted by 
~ 1180, ~ 790, and ~ 1080 km s~^ with respect to the narrow component in the spectra of 
MJM 3, 7, and 20. The Mg II hues in the spectra of the QSOs MJM 3, 7, and 20 have been 
decomposed into a narrow and a broad component. The broad component is found to be 
blueshifted by ~ 1180, ~ 790, and ~ 1080 km s~^ with respect to the narrow component in 
the spectra of MJM 3, 7, and 20. This suggests that either the red wings of the Mg II lines 
of MJM 3, 7, and 20 are affected by absorption features more strongly than the blue wings 
or the broad components of the Mg II lines originate in the outflowing winds. 

5 CONCLUSIONS 

We have derived soft X-ray spectral shapes and light curves of a nearly complete sample of 
22 QSOs. We also presented optical spectra of 12 QSOs from our sample. Our main results 
follows. 

(i) About 33% of the QSOs show a long term (~ 2 yr) soft X-ray variability while only one 
QSO MJM 10 shows rapid variability. 

(ii) Only three QSOs, MJM 7, 10, and 15 out of 22 QSOs (~ 7%) show indications for the 
presence of significant intrinsic absorption. The former two QSOs show excess absorption 
during both the observations of 1991 and 1993, while MJM 15 does not show excess absorp- 
tion during the observation 1993. 

(iii) The soft X-ray photon index of the QSOs in our sample ranges from 1.4 to 3.7. The 
average photon index of the sample is 2.40 ± 0.09 with a dispersion of 0.57. 
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(iv) The average photon index of the QSOs is found to flatten at higher redshift. This can 
be understood in terms of the redshift effect of mean intrinsic QSO spectra consisting two 
components - a soft excess component and a power-law. 

(v) Only one QSO MJM 10 out of 22 has been found to be an NLQSO. 

(vi) The broad component of the Mg II line in the spectra of MJM 3, 7, and 20 are found 
to be blueshifted by ~ 1180, ~ 790, and ~ 1080 km s~^, respectively. 
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